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Abstract 

We explored the motion of test particles near slowly rotating relativistic star having a uniform 
luminosity. In order to derive the test particle's equations of motion, we made use of the radiation 
stress-energy tensor first constructed by Miller and Lamb [3]. From the particle's trajectory ob- 
tained through the numerical integration of the equations of motion, it is found that for sufficiently 
high luminosity, "suspension orbit" exists, where the test particle hovers around at uniform angular 
velocity in the same direction as the star's spin. Interestingly, it turned out that the radial position 
of the "suspension orbit" was determined by the luminosity and the angular momentum of the star 
alone and was independent of the initial positions and the specific angular momentum of the par- 
ticle. Also found is that there exist not only the radiation drag but also "radiation counter-drag" 
which depends on the stellar radius and the angular momentum and it is this radiation counter- 
drag that makes the test particle in the "suspension orbit" to hover around at uniform angular 
velocity which is greater than that induced by the Lense-Thirring effect (i.e., general relativistic 
dragging of inertial frame). 

PACS numbers: 04.20.-q, 97.60.Jd, 95.30.Gv 
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I. INTRODUCTION 



Astrophysical accretion flow onto massive or compact stars is one of tfie major concerns 
in Astronomy and Astropliysics. In tlie present work, we explore the ground work for the full 
understanding of the accretion flow onto highly luminous slowly rotating relativistic stars. 
In the current treatment of the accretion process, the effect of radiation pressure on the 
inflow has not been fully addressed. Therefore in the present work, we attempt to include 
systematically the role of the radiation pressure in the accretion process. To this end, we 
explore the effect of the radiation pressure on the motion of a single test particle. Eventually 
we hope that this elementary study of ours will be extended to the case of actual accretion 
flow which can be thought of as consisting of large number of constituent single particles. 

We now begin with the summary of the present status of the research in the literature 
along this line. In their pioneering work, Abramowicz, Ellis, and Lanza jl| demonstrated 
that radiation from the highly-luminous non-rotating spherical massive star generates a 
"critical point" above the stellar surface, but they limited the motion of the particles to 
one-dimensional radial direction alone and did not deal with the case of rotating central stars 

n 

emitting isotropic radiation. Later on. Miller and Lamb |2| extended one- dimensional motion 
to two-dimensional one and pointed out that the trajectory of the particle is significantly 
affected by the radiation from the non-rotating star if the luminosity of the star is greater 
than ~ 1% of the Eddington luminosity. They, however, confined themselves to the spherical 
symmetric spacetime describing the non-rotation of the central star. Some time later, they 
considered [3| the effects of slow rotation of the central star and constructed the radiation 
stress-energy tensor describing the radiation field from the slowly rotating central star. They, 
however, failed to notice the emergence of critical radius (which corresponds to the critical 
point reported in Abramowicz et al.[l]) as they considered the set-up in which the luminosity 
is well below the Eddington's critical value. 

In the present work, in order to understand the effects of the radiation pressure on the 
accretion onto the highly luminous rotating relativistic stars in a rigorous and complete 
.anne, we employ the .adiat.on *e=.ene.gy tenso. wMch . given by MiUe. and Lamb 
as an elaboration on that given originally by Abramowicz et al. [JJ to describe the radiation 
emitted from the slowly rotating central stars, and derive the equation of motion. 

By integrating numerically the equations of motion, we realized that there exists the 
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"suspension orbit" (where the test particle hovers around the central star) that corresponds 
to the "critical point" in And it turns out that the radial position of this "suspension 
orbit" depends on the luminosity and the angular momentum of the central star alone, 
and does not depend on the initial position and the initial angular momentum of the test 
particles. 

In addition, it is realized that the test particle in the "suspension orbit" has uniform 
azimuthal velocity in the same direction as the star's spin motion. It is interesting to note 
that the uniform azimuthal velocity of the test particle at the "suspension orbit" is greater 
than that due to the Lense-Thirring effect. 

In section 2, we provide the derivation of the equations of motion from the known radia- 
tion stress-energy tensor and discuss the usual radiation drag terms and the newly discovered 
radiation counter-drag terms. In section 3, we present the results of the numerical integra- 
tions showing the emergence of the "suspension orbit" . In section 4, we analyze the motion 
of the particle in the "suspension orbit", and finally, in section 5, we end with the discussion 
of our results. 

II. EQUATIONS OF MOTION 

To derive the equations of motion that govern the trajectory of the test particles in the 
presence of radiation from the slowly rotating relativistic star, we assume that radiation 
source emits isotropically the radiation from the whole stellar surface and we shall employ 
the following metric obtained from the Kerr black hole metric ( 6|) by retaining its terms to 
the first order in Kerr parameter a to represent the spacetime exterior to the slowly rotating 
star, 

ds^ = g^i,dx^dx'' 

= -(1 - 2M/r)dt^ + (1 - 2M/ry^dr'^ + r^dO^ + r^sm'^e{d(f)'^ - 2ud(f)dt) (1) 

where oo = 2J/r^ is the Lense-Thirring angular velocity (I?!) arising due to the frame- 
dragging effect of a stationary axisymmetric spacetime, whch can be identified with the 



orbital angular velocity of the LNRF(Locally Non-Rotating Frame; see [8|]; {q]), M and J 
are the gravitational mass and angular momentum of the star respectively. We work in 
the geometric units where G = c = 1 (G is the gravitational constant and c is the speed 

n 

of light). Following Miller and Lamb |3[, we introduce a dimensionless angular momentum 



j = cJ/{GM^) and dimensionless velocity v =< v'^ > /c as a convenient measure of the 
rotation rate of the star and the rotation rate of the radiation source, respectively, where 
< f > is the appropriate average (see \^) of v'^ (azimuthal linear velocity) over the emitting 
surface visible from the test particle. The hat denotes the physical quantity measured in the 
LNRF. In this paper 'slow rotation' means j <C 1 and f <C 1 and we keep terms that are only 
first-order in j and v. For a neutron star with radius R 10 km, mass M f» 1.4Mq, and 
spin frequency Ug ^ 600 Hz(i.e., a millisecond pulsar), the dimensionless angular momentum 
j is approximately 0.2. 

The equations of motion (which actually is the geodesic equation) are given by 

a" = ^ (2) 

m 

where /° denotes a radiation force exerted by the radiation (or luminosity) on the test 
particle, m is the rest mass of the particle, and 

«° = ^ + ^''.uU'U^ (3) 

is the acceleration, with U" being the four- velocity of the particle, F^^^ = \g°'^{gi3^^u + gi3u,^ — 
giiv,p) being the Affine connection, and comma (,) denoting partial derivatives. 

For the sake of computational convenience, like in jsj, we assumed that the radiation 
scatters off the test particles and the momentum-transfer cross section a of the test particle 
is independent of energy (frequency) and direction of the radiation. Hence the radiation 
force /" due to scattering of the radiation is proportional to and in the direction of the 



radiation flux in the comoving frame (particle's rest frame), and is given by (see lOj) 



r = <yF^, (4) 

where is the quantity obtained by transforming the radiation energy flux T^^ measured 
in the comoving frame using the orthonormal tetrad e? associated with the particle's rest 
frame as follows, 

po -arpiO 

= -h'^pT^'^U,, (5) 

where /i^ = 5^ + U^Up is the projection tensor that projects onto each spacelike hyper- 
surface and T^'^ is the radiation stress-energy tensor in first-order (in j) Boyer-Lindquist 



coordinates (see j3|)(for j 7^ 0). 



According to [3|, the components of the radiation stress-energy tensor T"'^ as measured 
in LNRF are given by, 

T" ^ 27rJo(r)(l -cosao) (6) 
T**" ?5i 7r/o(r) sin^ ao 

T^4> ^ ^/o(r)J^(r)(cos=^ao-3cosao + 2) 

^ ^/o(r)(l-cos3ao) 
T'i ^ ll,(r)J{r) sin' ao 

T^^ ^ -/o(r)(cos2ao-3cosao + 2) 
3 

T^4> ~ ^/o(r)(cos^ao-3cosao + 2), 

where the subscript denotes the quantity for the case of non-rotating star, oq is an apparent 
viewing angle of the star seen by a locally static observer in Schwarzschild spacetime and 
is given by sinao = { 1-2 m/r ) ^ (^^^ El) radius of the star R > 3M, and Io{r) 

is the frequency-integrated specific intensity at the radial position r and is given by (see 
Appendix A in 

(1 - 2M/R) mM ( 

where the Eddington luminosity L'^^^ = AnmM/a is the luminosity of a spherically sym- 
metric source such that at infinity the outward radiation force balances the inward gravity 



see 



and L°° is the luminosity of the star as measured by an observer at infinity, and 



J^{r) is given by 

By transforming the above radiation stress-energy tensor T^'^ to the LNRF using tetrad 
e'i (which is given below) associated with the LNRF, the radiation stress-energy tensor T°'^ 
in the first-order Boyer-Lindquist coordinates (see p]) is obtained as, 

T"'^ = e^e^T^^ (9) 

where the tetrad associated with the LNRF are, 

e° = {l-2M/ry/^dt, 



= (1 - 2M/r)-i/2^r, (10) 
= rde, 

e = —2j — — sin 6 + r sin 6'(i0. 

We now focus our attention on the orbits confined to the equatorial plane {6 = ^,Ug = 0). 
The decomposition into each component of the equations of motion (5) in tensor form is 
given in the Appendix. It should also be noted that as the radiation stress-energy tensor 
(equation (6)) first constructed by Miller and Lamb [sl is valid up to j = 0.20, the equations 
of motion derived from it also valid within above ranges. 

Since background spacetime of equation (1) has a rotational isometry, Killing theorem 
states that there exists a rotational Killing field = 5^ such that the test particle's specific 
angular momentum / = g^j^rj'^U'' = g^^d'^U'^ = g^uU" = is conserved. Therefore, the 
azimuthal component (equation (23)) of the equations of motion in the Appendix governs 
the time evolution of the test particle's specific angular momentum and can be rewritten as. 



dr 



fir) 



3(1- 2M/r) 
L fjr) 
3(1- 2M/r) 



4j 



A{ao)U,' + ( 1 



2M 



(cosao sin^ ao)f^r 



A{ao)U^ + rJ{r)B{ao) + ( - ) J{r)B{ao)U^ 



Ut 



Lf{r) 
Lf{r) 



(2 sin^ aoWtU^ + 4j [ — \ (sin^ ao)U, 



J{r){sm'ao)Ur 



(11) 



where A(ao) = cos^ ao ~ 9cosao + 8, B(ao) = cos^ ao ~ 3cosao + 2, L = ( ] is the 
luminosity parameter, and /(r) = ^ (ilw/i^^ • 

We are now ready to envisage the features of the solution to the azimuthal component of 
the equation of motion given in equation (11). Obviously, however, the analytic solution to 
this coupled nonlinear ordinary differential equation (11) is not readily available. Therefore, 
we shall look for its numerical solution in the next section. However, even before that we can 
read off the essential features of the solution. To this end, we will interpret this azimuthal 
component of the equation of motion as the equation that determines the time evolution 
of the test particle's specific angular momentum as f/,^ = /. To summarize, this azimuthal 
component of equation of motion breaks into three parts: the first part is the line one of the 
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equation (11) and the two terms in this hne are responsible for the radiation drag, that is, 
the well-known Poynting Robertson effect since it is manifest that the overall sign of this 
line is negative definite and these two terms are linearly proportional to the test particle's 
specific angular momentum and the star's luminosity. The second part is the line two of 
the equation (11) and three terms in this line are responsible for the radiation counter-drag, 
that has the effect opposite to the Poynting Robertson effect since it is obvious that the 
overall sign of the line two is positive definite. Third part consists of line three and line four 
of the equation (11) and terms with opposite signs can be regarded as being responsible 
for, say, radiation drag for one sign and radiation counter-drag for the other. Now, in what 
follows, let us be more specific on the nature of each term in these lines. For instance, the 
first term in line two is due to the Lense-Thirring effect, that is, the dragging of inertial 
frame and the second and third terms in line two are due to the spin of the central star. 
Next, second term in line three is due to the Lense-Thirring effect and two terms in line four 
are due to the spin of the central star. 

It is interesting to note that although the radiation drag or the Poynting Robertson effect 
has been long-known, the nature of this effect has not been unveiled manifestly thus far. In 
equation (11) above, however, we were successful to quantify this Poynting Robertson effect 
by explicitly identifying the terms in line one which are responsible for the effect. In other 
words, by examining the equation of motion of a test particle placed in the background of 
luminous relativistic spinning star, we demonstrated manifestly that actually the Poynting 
Robertson effect takes place. Interestingly, however, this is not the end of the story. Re- 
markably, in the test particle's equation of motion, there are also terms which appear to 
be responsible for the effect just opposite to the Poynting Robertson effect, that is, terms 
in line two and some more terms in line three and four. This last point appears to imply 
that particularly for luminous relativistic spinning star, the new effect just opposite to the 
Poynting Robertson effect, which will be coined henceforth as the "radiation counter-drag", 
takes place as well. To the best of our knowledge, the counter-drag effect of this sort has 
never been reported in the literature so far. Therefore, in the following sections, we will solve 
the test particle's equation of motion numerically to construct and investigate quantitative 
solutions that will support our analysis of the features of the solutions stated above. What 
is more, based on both this numerical analysis and analytical approach, we will report on 
the emergence of the "suspension orbit" (which turns out to be the extension of the critical 
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point pointed out in earlier study [1] in the absence of star's rotation). 



III. NUMERICAL INTEGRATION 

In this subsection, we shall present the numerical solution to the equations of motion derived 
in section 2. We begin with the brief description of our treatment of this numerical analysis. 
Using the equations of motion derived in section 2, we have followed the trajectory of the 
particle in the presence of radiation from the slowly rotating star. We confine the motion 
of the particle onto the equatorial plane so that the polar angle component of the velocity 
Ue is set to be zero {Ug = 0). We assume that the star has uniform density, so the angular 
velocity of the star and the angular velocity of the LNRF at the stellar surface are given, 
respectively, by = |j and u = 2j Thus, the average azimuthal velocity v"^ of 

the radiation source as measured by an observer in the LNRF is calculated to be. 

Fig. 1 shows the trajectories of the particles starting at the position of r = 6M where 
the particles have three azimuthal velocities of 0.10(dotted curve), 0.25(dashed curve), and 
0.30(dash-dotted curve), respectively. Initial radial velocities of all the particles are equal 
to zero and the luminosity parameter of the star with R = 4M is L = 0.75. Solid line with 
radius of AM denotes the stellar surface, long dashed circle is the virtual circular orbit with 
radius of r = 6M (i.e., ISCO; Innermost Stable Circular Orbit). The starting point of all 
the particles are the same as (6M, 0) in Cartesian coordinates {x,y). The rotation of the 
star is counter-clockwise. Hence the particles start counter-clockwise rotational motion. 

As can be seen in Fig. 1, although the three particles start out in different azimuthal 
velocities, they end up being along the same circular orbit which we henceforth shall refer to 
the "suspension orbit" . It turns out that this "suspension orbit" lies in between the stellar 
surface(i.e., r = AM) and the ISCO (i.e., r = 6M). According to the numerical analysis, its 
radius is given by r = 4.66M. 

Fig. 2 shows the radial velocities of the particles in Fig. 1 inflowing toward the stellar 
surface from r = 6M as a function of the radius. Numerical integration demonstrates that 
the radial velocities of the particles at the "suspension orbit" vanish and the time rate of 
change of the radial velocities is also equal to zero, i.e., f/^ = and ^ = at the ' suspension 
orbit" . 
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x-axis(M) 



FIG. 1: shows the trajectories of the particles having three azimuthal velocities of 0.10(dotted 
curve), 0.25(dashed curve), and 0.30(dash-dotted curve), respectively when they initially co-rotate 
with the star (which is rotating counter-clockwise all the way) having angular momentum j = 0.1. 

Fig. 3 shows the azimuthal velocities of the particles in Fig. 1 as a function of the 
radius, and there we can figure out that the azimuthal velocities v'^ of the particles at the 
"suspension orbit" as measured by a locally static observer are all equal to v'^ ~ 0.056 in 
units of c(the speed of light) and the time rate of change of the azimuthal velocities is 
equal to zero, i.e., [/^ = constant and ^ = at the "suspension orbit". These indicate 
that the particles orbit at constant speed there. Remarkably, essentially the same is true 
for the case where the central star is not rotating. This point appears to indicate that 
the emergence of both the critical point for the case of non-rotating central star and the 
"suspension orbit" for the present case of rotating central star are indeed a generic feature 
that a high luminosity relativistic star exhibits. 

Thus far in Figs. 1 through 3, we have studied the case where the central star and 
the test particles "co- rotate" all the way. Next, we move on to the other case where they 
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FIG. 2: shows locally measured radial velocities as functions of radius for particles in Fig. 1 
inflowing from r = 6M toward the "suspension orbit" . 

counter-rotate initially but end up co-rotating eventually and here, the rotation of the star 
is counter-clockwise all the way. The result is given in Figs. 4 through 6. 

To summarize, in Figs. 1 through 3, we have studied the case where the central star and 
the test particles "co-rotate" all the way. In Figs. 4 through 6, however, we have studied the 
other case where they counter-rotate initially but end up co-rotating eventually. Namely, 
regardless of initial conditions, i.e., whether they are initially co-rotating or counter-rotating, 
the system reaches the same final equilibrium state where the test particles end up co- 
rotating with the central star. This result indeed is very interesting and curious particularly 
for the case when the central luminous star is spinning since the eventual fate of the test 
particles is the co-rotation with the central star at the "suspension orbit". Therefore, we 
need careful understanding of the underlying physics and our interpretation is based upon 
the geodesic equation of the test particles at the "suspension orbit" given in equation (13) 
below and it can be describes as follows. 

Firstly, first line term in equation (13) plays the role of radiation drag because its overall 
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FIG. 3: shows locally measured azimuthal velocities as functions of radius for particles in Fig. 1 
inflowing from r = 6M toward the "suspension orbit" . 

sign is negative definite due to Ut < and is proportional to the test particle's specific 
angular momentum U^. Note that the radiation drag of the first line term works on the test 
particle regardless of the spin of the central star. 

Secondly, fourth line terms, on the other hand, play the role of "radiation counter-drag" 
because its overall sign is positive definite and is proportional to the central star's angular 
momentum j. Note that the radiation counter-drag of these fourth line terms work on the 
test particle regardless of the test particle's specific angular momentum U^. 

Indeed, the radiation drag term in the first line of equation (13) that represent the 
Poynting-Robertson effect is well-known and it is generic as it is independent of the central 
star's angular momentum. The radiation counter-drag terms in the fourth line has been 
neglected in our conventional understanding of the Poynting-Robertson effect and hence is 
rather unfamiliar. Besides, it is not so generic as it appears only in the presence of the 
central star's angular momentum. Therefore, unlike the case where the central luminous 
star is non-rotating, for the case at hand where the central star is spinning, not only the 
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FIG. 4: shows the trajectories of the particles having three azimuthal velocities of -0.10(dotted 
curve), -0.25(dashed curve), and -0.30(dash-dotted curve), respectively when they initially counter- 
rotate with the star(which is rotating counter-clockwise all the way) having angular momentum 
j = 0.1 but end up being co-rotating with the star mainly due to the radiation counter-drag. 

radiation drag term in the first line but the radiation counter-drag terms in the fourth 
line operate as well and indeed they would affect the motion of the test particle on equal 
footing to determine its trajectory. Interestingly enough, our numerical study above exhibits 
that once the test particle arrives at the "suspension orbit" , it never comes to a complete 
stop. Rather, it keeps co-rotating with the central luminous spinning star and this is the 
new discovery of the present work which has not been realized in the previous literature 
addressing the similar related issues. 

As can be noticed from the numerical integration, the radiation from the slowly rotating 
star makes the test particles to hover around the star with uniform azimuthal velocity 
regardless of the initial position and the initial angular momentum of the particle, and the 
particle's motion in the "suspension orbit" is characterized by the following conditions, 

Ur ^ 
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FIG. 5: shows locally measured radial velocities as functions of radius for particles in Fig. 4 
inflowing from r = 6M toward the "suspension orbit" . 



dUr 



constant 



dr dr 

IV. EXAMINATION OF THE MOTION IN THE "SUSPENSION ORBIT" 

In order to understand the nature of "forces" exerted on the particle hovering around the 
"suspension orbit", let us examine the equations of motion using the conditions mentioned 
above in section 3 and then determine the coordinate radius of the "suspension orbit" and 
the azimuthal velocity of the particle as measured by the locally static observer. 

By inserting Ur — and ^ = of the conditions for the "suspension orbit" into the 
0-component (equation (11)) of the equations of motion, we obtain the following. 



M2 

4j—A{ao) + 2J{r)B{ao) 



Ul + rA{ao)U^Ut + r^J{r)B{a^) = 0, 



(12) 



where the luminosity parameter L as a common factor is omitted. Then, plugging J{r) in 
equation (8) into the above equation (12) gives, after some manipulation, 
dU^ 



= 



dT 



= A{aQ)UtU^ 
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FIG. 6: shows locally measured azimuthal velocities as functions of radius for particles in Fig. 4 
inflowing from r = 6M toward the "suspension orbit" . 



Firstly, terms in line one in equation (13) play the role of radiation drag because its 
overall sign is negative definite due to C/^ < and is proportional to the test particle's specific 
angular momentum C/^. Therefore, the radiation drag terms in line one work on the test 
particle having the azimuthal velocity regardless of the spin of the central star. Secondly, 
as [A(ao) — 4i?(Q;o)] of hne two in equation (13) has positive value for ^ 0, the overall 
sign of line two is positive definite, thus this line serves as the radiation counter-drag that 
speeds up the azimuthal motion. It should also be noted that this line two is proportional 
to the Lense-Thirring angular velocity u — 2j (^^^ and the test particle's specific angular 
momentum square C/|, thus the radiation counter-drag of this libe two is due to the Lense- 
Thirring effect arising from the rotation of the central star. Unfortunately, this line two in 




(13) 
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equation (13) has very tiny contribution to the azimuthal velocity of the particle because of 
f/(^ ~ j <C 1 in the "suspension orbit" . This proportionality relation between and j can 
be noticed from the fact that the "suspension orbit" in slow rotation case amounts to critical 
point (see [l|) in non-rotating case and at the critical point is equal to zero, and also can 
be supported by subsequent calculation in equation (15). Therefore, all terms including 
is negligible in the "suspension orbit" . However, in the case where it is not allowed to ignore 

like in initial trajectory at r = 6M, all the radiation counter-drag terms including f/| 
could contribute somewhat to the azimuthal velocity of the particle. Thirdly, two terms in 
line three in equation (13) have overall positive definite sign, and thus line three behaves as 
the radiation counter-drag. Since this line three is proportional to the cental star's angular 
momentum j, the radiation counter-drag of line three is due to the rotation of the central 
star. Also, as line three is proportional to f/|, its contribution is negligible in the "suspension 
orbit" . Lastly, two terms of line four in equation (13) have overall positive definite sign, thus 
line four serves as the radiation counter-drag. Since line four is linearly proportional to the 
cental star's angular momentum j like the third line, the radiation counter-drag of line three 
is also due to the rotation of the central star. However, the radiation counter-drag of the 
fourth line acts on the test particle regardless of the velocity of the particle. Also, they are 
comparable to the radiation drag term of line one. Therefore, the uniform azimuthal velocity 
of the particle in the "suspension orbit" can be attributed to the terms of line four. In other 
words, the balance between first line term (radiation drag) and fourth line terms (radiation 
counter-drag) makes the particle to hover around the central star with uniform azimuthal 
velocity. Since the radiation drag term of the first line exerts on the test particle having the 
specific angular momentum f/^ regardless of the spin (j) of the central star, this term exists 
even when the central star is not rotating, whereas the radiation counter-drag terms in line 
two through four are linearly proportional to the central star's angular momentum (j), if 
the central star is non-rotating, all the radiation counter-drag terms disappear. Therefore, 
the emergence of the radiation counter-drag is attributed to the rotation of the central star. 
It also is of interest that there is another radiation counter-drag term proportional to the 
Lense-Thirring angular velocity u and the particle's specific angular momentum square (f/|). 

Next, in order to calculate particularly the azimuthal velocity as measured by a locally 
static observer, we make use of normalization condition g^'^U^Up = — 1. That is, plugging 
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Ur — oi the conditions for the "suspension orbit" into the normahzation condition yields, 

(l -—)Ul- Aj—U,U, + r^(l-—y r^Vl = 0. (14) 

It is difficult to obtain analytically f/^ and \Jt from the combination of equations (12) 
and (14). Therefore, by assuming that [/| is negligible from our experience in the case of 
non-rotating star, i.e., C/| ~ and neglecting terms of order higher than linear in j', we can 
get approximate and respectively, as. 



[/, ^ r (l - ^) J{r) 



A{a) 



(15) 



— .P)-.-(l-^r- (16) 

Using i7(r) ~ j in equation (8), equation (15) indicates C/^ ~ j which is consistent with 
our assumption above. 

The azimuthal velocity of the test particle as measured by a locally static observer is 
given by 

1 U-e^ 



U -e 



where e'^ and e* are tetrad associated with the locally static observer. Substitution of 

^ ^ -r(l- Jir)^^ obtained from equations (15) and (16) into equation (17) 



Ut \ r ) ^ \ ' A{ao) 

yields the azimuthal velocity of the particle in the "suspension orbit" as measured by the 
locally static observer to be. 



^au-yl-— W^ + 2j — ) ■ (18) 



First term in equation (18) indicates the contribution to the azimuthal velocity of the 
particle in the "suspension orbit" due to the combination of the rotation of the star and the 
rotation of the radiation source which rotates with the star by being attached to the stellar 
surface. Since second term in equation (18) is equal to the azimuthal velocity of the LNRF 
measured by the locally static observer, it is the contribution due to the Lense-Thirring 
effect (dragging of inertial frame) arising from the slow rotation of the star. 
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In order to obtain the coordinate radius from the center of the star to the "suspension 
orbit" , we make use of ^ = and Ur — oi the conditions characterizing the "suspension 
orbit" together with C/^ in equation (15) and Ut in equation (16). By inserting these into 
equation (22) in the Appendix, that is the r-component of the equations of motion, we get 



r OO 



2M\^/2 



and the solution of equation (19) is the coordinate radius r suspension from the center of 
the star to the "suspension orbit". As we can notice from equation (19), it is noteworthy 
that the coordinate radius Vgo of the "suspension orbit" indicating the circular orbit of 
the particle hovering around the slowly rotating massive star depends on the luminosity 
parameter ( ) and the dimensionless angular momentum j of the star alone, and does not 

\^Edd J 

depend on the initial position and the initial velocity (angular momentum) of the particles. 
Also, since equation (21) in the Appendix , that is the i-component of the equations of 
motion, is satisfied by equations (12) and (14), the problem of over-determining Ut and 
does not occur. For ( ) = 0.75 and j — 0.10, equation (19) gives Vgo — 4.63M which is 

\^Edd J 

similar to the value {vso — 4.66M) in the numerical integration of section 3. Then, plugging 
Vso ^ 4.63M and j = 0.10 into equation (18) yields fftispension — 0.056 = 0.044 + 0.012 which 
is also very similar to the value {v'^ ~ 0.056) obtained from the numerical integration in 
section 3, where 0.044 is the contribution from Doppler shift resulting from the combination 
of the rotation of the star and the radiation source, and 0.012 is the contribution from the 
Lense-Thirring effect due to the rotation of the star. Thus, we can find that the Doppler 
shift due to the rotation of both the central star and the radiation source have even larger 
contribution to the azimuthal velocity of the test particle than Lense-Thirring effect by the 
rotation of the central star. 

In Fig. 7, as can be seen from the comparison of the solid line (the plot of the approximate 
expression given by equation (19)) and numerical values (triangles and asterisks) obtained 
by the numerical integration for j — 0.10, above approximate equation (19) is valid only for 
j < 0.10. li j — and v — 0, i.e., the star is non-rotating and the radiation source does not 
rotate either, equation (19) reduces to, 

f \ ( 2M\^/^ 



r OO 



(l--) . (20) 
In this case, for (:^) = 0.75, equation (22) gives Tso — "i^M/l ~ 4.57M which is 
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FIG. 7: shows the coordinate radius of the "suspension orbit" rgo as a function of the luminosity L 
for the angular momentum of the central star j = 0.10 and 0.20, respectively. Solid line denotes the 
plot of the approximate expression given by equation (19) for j = 0.10, and Triangles (j = 0.10) 
and asterisks (j = 0.20), respectively, denote values obtained through the numerical integration. 

smaller than Vgo — 4.66M in the slowly rotation case. Thus, the rotation of the star and 
the radiation source makes the "suspension orbit" to expand outward, which is obvious 
from the fact that the rotation of the particles generates centrifugal force and this force is 
balanced with the gravitational force and the radiation pressure. 

Also, if Vso = 4Af (stellar radius) and Vgo = QM (ISCO) in the case of j = 0.10, then 
from equation (21) we can get ( ] — 0.70 and ( -f^ ] — 0.80, respectively, so if the 

\^EddJ \^Edd J 

luminosity of the star measured at infinity hes within the range of 0.70 < ( ] < 0.80, 
the "suspension orbit" appears in boundary layer between the stellar surface and the ISCO. 
Therefore, if the slowly rotating central star with radius of i? = 4M has luminosity of 
L = 0.75 and angular momentum of j = 0.10, the test particle hovers around the star with 
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azimuthal velocity of about 0.05c in the boundary layer between the stellar surface and the 
ISCO where the rotation direction of the particle is the same as that of the star's spin. 



V. DISCUSSION 

We now summarize what we have realized in the present work which are new ingredients 
that have not been addressed in the previous literature. The results presented in sections 
3 and 4 show that there exists "suspension orbit" whose location is determined by the 
central star's luminosity L as measured at infinity and the star's angular momentum j and 
is independent of the initial positions and velocities of test particles. To be more concrete, 
we have explored the 2-dimensional motion of test particles on the equatorial plane around 
slowly rotating star with the asymptotic luminosity measured at infinity in the range of 
0.70 < ( ) < 1- We found out that for the luminosity in the above range, there exist 

\^Edd J 

"suspension orbit" at which the radial velocity of the test particle vanishes and proper time 
rate of change in the radial and azimuthal velocity also vanish and the particle's azimuthal 
velocity is constant, thus the test particle hovers around the star with uniform azimuthal 
velocity regardless of their initial positions and velocities. And it is interesting to note 
that not only the radiation drag but also the radiation counter-drags which result from the 
central star's spin exert on the particle in the "suspension orbit" and the balance between 
the radiation drag and the radiation counter-drags makes the particle to hover around the 
star at uniform azimuthal velocity much greater than that due to the Lense-Thirring effect 
(i.e., the dragging of inertial frame). Furthermore, it is noticeable that there exists another 



radiation counter-drag term which is proportional to t 



le Lense-Thirring angular velocity u. 



Interestingly enough, Bini, Jantzen, and Stella Uj reported a study which happens to 
be akin to the motivation of our work and consistent with our results in some respects. 
They found the existence of the critical radius at which the test particle co-rotates with the 
geometry in the Kerr spacetime background. However, the radiation stress-energy tensor 
they employed is valid for only the photons in outward radial motion with zero angular 
momentum, thus is not appropriate to the applications to the accretion process onto the 
rotating relativistic stars like the neutron star. In the present study, on the other hand, 
since we employed the radiation stress-energy tensor first constructed by Miller and Lamb 
[sl which has no limitations whatsoever on the character of the emitted photons, our results 
can be applied to the accretion process onto the rotating relativistic stars. 
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Lastly, in the forthcoming article, we will report on our result of the study where the 
central luminous star is non-rotating {j — 0). 
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Appendix 



The equations of motion (5) of a tensorial form is decomposed into each components as 
follows: 



dr 
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where A{ao) = cos'^ ao — Ocoscto + 8, B{ao) = cos^ cto — 3cosq;o + 2, and sincto 
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